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Abstract 



Tfiis paper is a detailed description of an algorithm based on a gener- 
alized Buchberger algorithm for constructing Grobner-type bases asso- 
ciated with polynomials of shift operators. The algorithm is used for 
calculating Feynman integrals and has proven itself efficient in several 
complicated cases. 
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1 Introduction 

Feynman integrals over loop momenta are building blocks of quantum-theoretical 
amplitudes in the framework of perturbation theory. After a tensor reduction based 
on some projectors a given Feynman graph generates various scalar Feynman integrals 
that have the same structure of the integrand with various distributions of powers of 
propagators which we shall also call indices: 



Here ki, i = 1, . . . ,h, are loop momenta and the denominators Er are either quadratic 
or linear with respect to the loop momenta Pi = ki, i = 1, . . . ,h, or independent 
external momenta ph+i = qi, . . . ,Ph+N = In of the graph. (Linear propagators 
appear naturally in effective field theories and in asymptotic expansions of Feynman 
integrals in various limits.) 

Irreducible polynomials in the numerator can be represented as denominators 
raised to negative powers. Usual prescriptions k"^ = k'^ + iO, etc. are implied. The 
dimensional regularization pQ with ci = 4 — 2^ is assumed. The Feynman integrals 
are functions of the masses, m^, and kinematic invariants, qi ■ qj. However, we shall 
omit this dependence because we shall pay special attention to the dependence on 
the indices. We shall also omit the dependence on d. 

At the modern level of perturbative calculations, when one needs to evaluate 
millions of Feynman integrals (^, a well-known optimal strategy here is to derive, 
without calculation, and then apply some relations between the given family of Feyn- 
man integrals as recurrence relations. A well-known standard way to obtain such 
relations is provided by the method of integration by parts (IBP) j2]. Practically, one 
starts from IBP relations 



After the differentiation, resulting scalar products, fcj ■ kj and ki ■ qj are expressed in 
terms of the factors in the denominator, and one arrives at IBP relations which can 
be written as 



where bij are integer, Cj are polynomials in aj, d, masses rrii and kinematic invariants, 
and F{ai, . . . , a„) are Feynman integrals (0) of the given family. 

One tries to use IBP relations in order to express a general dimensionally regular- 
ized integral of the given family as a linear combination of some 'irreducible' integrals 
which are also called master integrals. There are several recent attempts to make 
this reduction procedure systematic, in particular, the so-called Laporta's algorithm 
H] (There is a public version of implementing the corresponding algorithm on a 








(3) 
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computer and Baikov's method (see jHllZl and chapter 6 of |S]). Another activity 
in this direction is connected with the use of Grobner bases [U^. The first variant of 
this approach was suggested in ^Hj, where IBP relations were reduced to differential 
equations. First attempts to use directly the non-commutative Grobner bases in the 
algebra generated by shift operators were made in [HI [T2] . In the previous paper [13] 
(see also for a brief review) we presented another approach based on Grobner 
bases. Its specifics lies in using more information on the given family of Feynman 
integrals, in particular the boundary conditions, i.e. the conditions of the following 
form: 

F(ai, 02, ... , an) = when < 0, . . . a^^ < (4) 

for some set of indices ij (for example, we always have F(ai,a2, . . . , a„) = if all 
Oj are non-positive). Since that time the algorithm has been made more efficient by 
introducing the so-called s-form. Moreover, the algorithm has been made about ten 
times faster. Here it will be described in details (the paper [131 designed to give 
an introduction to the method and not all definitions were formal). 

2 Preliminaries 

To describe the algorithm we have to introduce some notation. Let /C be the field of 
rational functions of physical variables rrii, qi ■ qj, d, and A be the algebra^ over /C 
generated by elements Yi, Y[~ and Ai with the following relations: 

YrYr = YfYr , YrY, = Y^r^ Yr A, = A^r - 5,,^,, YrY., = l 

where = 1 if i = j and otherwise. For convenience we will write (yj )^ = Y~^ . 
Let T be the field of functions of n integer arguments oi, 02, . . . , a„. The algebra A 
acts on this field'^, where 

{Yi ■ F){ai, a2, . . . , a^) = F{ai, aj_i, + 1, a^+i, . . . , a„) , (6) 
{Ai ■ F){ai, 02, ... , an) = aiF{ai, 02, ... , a„) . 

The left-hand sides of relations ^ can be represented as elements of the algebra 
A applied to F; we will denote these elements by /i, . . . , fk- Now, for F(ai, . . . , an) 
defined by (H)), we have 

/,-F = Oor (/,-F)(ai,...,a„) = (7) 

^An algebra over a field is a vector space over this field and a ring at the same time. 

^{i) for any a € A and / G we have an element a ■ f € !F; {ii) for any a,h & A and f,g^J- 
we have (a + b) ■ {f + g) — a-f + a- g + b- f + b-g; (iii) for any a,b £ A and f E J- we have 
(a6)-/ = a- (&•/). 
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for all i. Let us generate a (left) ideal I by the elements fi, . . . , fk. We will call X 
the ideal of the IBP relations. Obviously, 

/■F = 0, or (/■F)(ai,...,a„) = for any fel. (8) 

Let us consider the set V with elements {ci, C2, . . . , c„} where all Cj are equal 
to 1 or —1. The elements of this set will be called directions. For any direction 
u = {ci, . . . , c„} we will consider a region = {(ai, • • • , an) ■ (a^ — l/2)cj > 0} and 
call it a sector. Obviously the union of all sectors contains all the integer points in 
the n-dimensional vector space and the intersection of any two sectors is an empty 
set. For a sector cr^ we will say that its direction is u. 

We will say that an element X G ^ is written in the proper form if it is represented 

as 

X = J2rMu---,An)llYf^'\ (9) 

i 

where vj are polynomials (with coefficients in /C) and dij are integers. (So, all the 
operators Ai are placed on the left from the operators Yi.) Obviously any element 
X & A has a unique proper form. We will say that an element of ^ is a monomial if 
in its proper form only one coefficient function rj is non-zero. 

Let = {(61, . . . , bn)} where all fej are integers and bi > 0. This is a semi-group 
(with respect to (61, ... , bn) + {b[, . . . , b'^) = {bi + b[, . . . , bn + b'^)) . We will say that 
an ordering on is proper if 

i) for any a G not equal to (0, . . . 0) one has a y {0, . . .0) 

ii) for any a,b,c & one has a y b if and only if a + c >~ b + c. 

Let us fix a direction z/ = {ci, C2, . . . , c„} We will say that the v-degree of a 
monomial r{Ai, . . . , An) Hi is {diCi, . . . , dnCn} if all the products diCi are non- 
negative and undefined otherwise. 

Let us choose a proper ordering. Take an element X & A and write it in its 
proper form. Consider the set of i/-degrees of all the monomials in the decomposition 
of X. If this set is empty (all the degrees are undefined) we will say that the i/-highest 
degree of X is undefined as well, otherwise we will say that i/-highest degree of X is 
the highest element of that set in terms of the fixed ordering. 

We will say that a direction {ci, . . . c„} is lower than {c'^, . . . c^} if ci < c'^, . . . , c„ < 
and they are not equal. The same is said about the corresponding sectors. We will 
say that a sector a is trivial if all the integrals -F(ai, . . . , a„) are zero for (ai, . . . , a„) G 
a due to boundary conditions (jH). The same will be said about the direction of the 
sector. 

3 The s-reduction algorithm 

Suppose that for each non-trivial sector a^, we are given a finite basis {X^ i, . . . , X^^^^} C 
X and an ordering. 
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We are going to describe an algorithm called s-reduction. 
Input: a linear combination of integrals . . . , a^). 

Output: another linear combination of integrals containing the integrals that 
could not be reduced by these bases (an integral is called F(ai, . . . , On) irreducible 
for the given sets of orderings and bases if the z/-reduction of F{ai, . . . , a^) returns 
F{ai, . . . , a„), where u is such a direction that (ai, . . . , a„) G a,y). 

If one chooses properly the orderings and the elements X then there is a finite 
number of irreducible integrals, therefore the s-reduction will always output a com- 
bination of them but such a choice is a non-trivial procedure. The algorithm that 
might construct such elements will be described in the next section. But first we 
have to explain what s-reduction is. Some examples of the bases can be found at 
http://www.srcc.msu.ru/nivc/ ab out /lab /lab4_2 / index_e ng . htm| 

S-Reduction 

1. L=Input; M = 0. 

2. While L^O 

3. Let 5* be the set of sectors that contain a point (ai, . . . , a„) where 
F{ai, . . . , a„) has a non-zero coefficient in the decomposition of L. 

4. Let 1/ be a direction such that E S and there is no other sector 
(7^1 G S such that ay is lower than a^i. 

5. Let L = Li + L2 where Li contains those and only those 
F(ai, . . . , a„) where (ai, . . . , a„) G a^. 

6. N = i^-reduction of Li by {X^^i, . . . , Xy^k„} (to be described below) 

7. Let N = Ni + N2 where Ni contains those and only those 
F{ai, . . .,an) where (ai, . . . , a„) G a^. 

8. L = L2 + N2; M = M + Ni. 
9. Output=M 

In a few words we are using the zz-reduction of X starting from higher sectors. 
Basically, it is a formalization of the standard method people use to solve the IBP 
relations "by hand" . Now we are going to describe what zz-reduction is. 

u-reduction 

Input: direction u; an element L that is a linear combination of integrals F{ai, . . . , an) 
such that all (ai, . . . , a„) G 0",^ ; a finite basis {Xi, . . . , X^} C X. 

Output: another linear combination of integrals F(ai, . . . , a„) such that (ai, . . . , a„) G 
a^/ where u' = u or v' is lower than v. 

1. Let p = {pi, ...,pn} = {(ci + l)/2, . . . , (c„ + l)/2}, where u = {ci, . . . ,Cn}. 

2. Set Y = 0. 

3. Let X G ^ be the element obtained by replacing F(ai, . . . , a^) 
with Hi Y^'^'-P' in L. One has L = {X ■ F){p). 

4. WhileX^O 

5. Let U be the highest term oi X,U = C UiYf'''\ 
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6. For all possible products T = {UiYj^') ■ Xj such that the proper form of T 
has a non-zero coefficient at UiYf''^' 
Do (there is a finite number of possibilities) 

7. Replace all Ai in the proper form of T with pi. 

8. If for some j such that Cj = —1 there is a term of T with Yj 

in a positive degree then continue the cycle with the next element. 

9. If the coefficient C of T at UiYf''^' is equal to zero 
then continue the cycle with the next element. 

10. Take Z = X — {C/C')T. Let Zi be the //-sector part of Z 
(the sum of monomials in the decomposition of Z 

that have a defined degree) and Z2 = Z — Zi. 

11. If the i/-highest degree of Zi is lower than the 
//-highest degree of X then replace X with Zi , 

Y with Y + Z2 and go to the start of the While cycle (step 4). 
12. Replace X with X -U andY with Y + U 
13. Return (Y ■ F){p) 

The idea of the algorithm is to represent the given linear combination as an element 
X of A being applied to F and the value being taken in the corner of the sector. Then 
one tries all possible transformation of the obtained element that lower the //-degree 
of X (this is a generalization of the standard reduction procedure). The point is that 
through all the algorithm the value {{X + Y) ■ F){p) is not changed. This is based on 
the fact that {Z ■ F){p) = for any Z e X. Note that the elements X and Y in their 
proper forms do not contain the operators Ai (because of the replacement in step 7). 
Therefore in step 9 we have C" G /C and C 7^ 0, so the division is possible. The step 
12 is the place where the terms that cannot be reduced are moved from X to Y. If 
the basis is chosen properly, those terms will correspond to a finite number of master 
integrals. 

4 Constructing s-bases 

Let us fix a non-trivial direction z/ = {ci, C2, . . . , Cn} and let p = {pi, . . . ,pn} = 
{(ci + l)/2, . . . , (c„ + l)/2}. Our task is to construct an s-hasis for this direction, i.e. 
such a basis {Xi, . . . ,Xfc} C 1 that the //-reduction for this sector has only a finite 
number of irreducible integrals. 

The relations Q give us a basis of 1 but generally it is not an s-basis. Moreover, 
one has to choose an appropriate ordering. 

We are going to describe an algorithm that takes the relations (jS)) as input and 
aims to construct an s-basis. The results of the algorithm greatly depend on the 
choice of the ordering and in complicated cases the proper choice seems to be unique. 
Now suppose we have fixed an ordering and an initial basis {Xi, . . . C X. Let 
us describe our algorithm. 
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First of all let us define the s-form of an element X e A (note that this definition 

and some more definitions below depend on the choice of the direction u, but it is 
fixed in this section). So, the s-form of an element X e ^ is an element T of the 
form (JliY^^^) ■ X satisfying the following properties: 

(i) The i/-highest degree of T is defined and for any integer {yi, . . . ,yn) such that 
ViCi > 0, . . . , ynCn > the i/-highest degree of (UiY^^) • T is equal to the i/-highest 
degree of T plus (t/i, . . . , 

(ii) The //-highest coefficient of (IliY^')-T does not vanish when Ai — pi, . . . , A^ — Pn 
for all integer (yi, . . . , y^) such that yiCi > 0, . . . , y„c„ > 0; 

{in) For all j such that cj — —1 the degrees of Yi in the proper form of T are non- 
positive; 

{iv) The numbers {ciXi, . . . , CnXn) are minimal possible for all {xi, . . . , x„) satisfying 
the properties {i)-{m). 

Let us reformulate this definition less formally. The properties {i), {ii) mean that 
this element has "enough" terms whose degrees lie in the sector cr^. It is needed so 
that this element can be used for the i/-reduction. The property (iii) means that 
if Cj = 1 then both Yj and Yj~ can appear in the s-form, but if Cj = —1 then only 
Y~ can. The meaning of this requirement is to control that using this element for 
//-reduction we will not obtain elements that lie in sectors higher than v. 

Through all the algorithm we will store a basis of X consisting of elements in the s- 
form. Let us describe how the reduction of an element of A via a basis {Xi, . . . , X^} C 
I works. 

1. y= s-form of Input. 

2. If y = then Return 0. 

3. For all j such that all the numbers di — d[ are positive, where (di, . . . , d„) is 
the //-highest degree of Y and {d[, . . . , d^) is the i/-highest degree of Xj 
Do 

4. Let Z = C'Y -C niF/'^'^'"'^^^ ■ Xj, where C is the //-highest coefficient 

of Y and C is the //-highest coefficient of (IljF-'^'^'^' '^'^) • Xj. 

5. Z = s-form of Z. 

6. If the i/- highest degree of Z is lower than i/-highest degree of Y then 
replace Y with Z and Go to step 2. 

8. Return Y. 

It is easy to see that this reduction stops after a finite number of steps (there 
can't be an infinite sequence of decreasing degrees). Basically, this procedure is close 
to the standard reduction procedure in the Buchberger algorithm. The difference is 
the usage of s-forms and the fact that the elements Yj can have both positive and 
negative degrees. 

Now we can describe the main algorithm. As it has been said earlier, it starts 
from a basis of X, moreover, all elements are taken in their s-forms. The goal of the 
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algorithm is to construct another basis such that the i/-reduction for this sector has 
only a finite number of unreducible integrals. Therefore, after an clement is added 
to the basis or a basis element is replaced we are performing a test to verify this 
condition (Completion Criteria). It consists of checking, whether for any j there is 
such m that for all integer / > m the element F{pi, . . . ,pj_i,pj + lcj,pj-^i, . . . ,Pn) 
can be reduced via this basis in the i/-reduction algorithm. 
Here is the main algorithm: 

1. While not Completion Criteria 

2. If there is an element in the basis that can be reduced by some other 
clement, replace it with the result of the reduction and restart the cycle. 

3. Choose a pair of elements of the basis X' and X". 

4. Choose the smallest possible integers {di, . . . , dn) such that dj > d'j and 
dj > d", where (li'^, . . . , d^) is the i/-highest degree of X' and (d'/, . . . , d^) 
is the i/-highest degree of X". 

5. Evaluate Z = C" UiY^'^'^*"^'^ ■ X' - C UiY^'^'^*~'^'^ ■ X", where C is the 
i/-highest coefficient of (IljF/*^* '^'^) • X' and C" is the i/-highest 

coefficient of (HiF/^^*"'^*'^) • X", — the s-polynomial of X' and X". 

6. Z = Reduction of Z. 

7. If Z 7^ then add Z to the basis. 

When implementing the present algorithm it is natural to store the information 
about the pairs where the s-polynomials have been evaluated to avoid repeating the 
same calculations. Of course, the choice performed at line 3 might be different, and 
the algorithm effectiveness greatly depends on this choice. One more improvement 
of the algorithm is the use of symmetries of the diagram. Instead of evaluating an 
s-polynomial one might take an element symmetric to some element of the basis and 
reduce it the same way. All this is realized by introducing a function on the set of 
pairs of elements of the basis and another one on the set of all possible symmetric 
element to the elements of the basis. So, at line 3 we are choosing such an element to 
evaluate, that the value of this choice function is minimal. Currently the algorithm 
uses the choice function intended to minimize the degree of the resulting element, 
and it already makes the algorithm effective. The work on finding the optimal choice 
functions is in progress. 

5 Conclusion 

The algorithm described above is close to the Buchberger algorithm, but has a signif- 
icant difference. The main problem is that generally one can have two elements of the 
same degree that cannot be reduced one by another in this algorithm. Therefore, the 
standard proof that the algorithm has to stop at some point does not work here. To 
prove that this algorithm stops for any given family of Feynman integrals is an open 
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problem which, hopefully, will be solved in the nearest future. On the other hand, the 
present algorithm appears to be much more efficient to solve reduction problems for 
Feynman integrals than more or less straightforward generalization of the Buchberger 
algorithm (see, e.g., This was demonstrated in [13j where the reduction prob- 

lem was solved for Feynman integrals relevant to the two- and three-loop static quark 
potential, with the number of indices n = 7. (In two loops, well-known results [T3] 
were reproduced.) New results obtained with the help of this algorithm, in cases with 
the number of indices n = 9, will be published soon jTHI. Preliminary analysis shows 
that the algorithm can work successfully in problems with the number of indices up 
to n = 12. The extension to more complicated problems will, of course, require its 
modifications, in various places. 

An implementation of the s-reduction part of the algorithm in Mathematica can 
be found at http:/ /www. srcc.msu.ru/nivc/about/lab/lab4_2/index_eng. htm together 
with pedagogical examples and at least the s-bases that have been constructed for 
the problems considered in |113'j . 
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